The ground state entanglement in the XXZ model 
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In this paper, we investigate spin entanglement in the XXZ model defined on a d-dimensional 
bipartite lattice. The concurrence, a measure of the entanglement between two spins, is analyzed. 
We prove rigorously that the ground state concurrence reaches maximum at the isotropic point. For 
dimensionality d > 2, the concurrence develops a cusp at the isotropic point and we attribute it to 
the existence of magnetic long-range order. 
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Entanglement, as the exhibition of pure quantum cor- 
relations between separate systems, has become one of 
the trademarks of the quantum mechanics for its non- 
local connotations 0, 01 ■ Recently, many physicists have 
made great efforts to understand the quantum entangle- 
ment in the ground states of some many-body spin mod- 
els [3lll.l5lla[7l.8...9....1Q..ll.]. One expects that a thorough 
investigation on the entanglement in these systems will 
provide new insight into the quantum phase transition in 
these systems |,12j . For example, Osterloh et a/j3] stud- 
ied the concurrence, a measure of entanglement of two 
qubitsfl^. between two spins located on a pair of nearest- 
neighbor sites in the transverse-field Ising model 4] . They 
found that this quantity shows singularity and obeys the 
scaling law in the vicinity of the quantum phase transi- 
tion point of the system. On the other hand, for other 
models, such as the antiferromagnetic XXZ chain, the 
concurrence behaves in a completely different wayj2j. As 
shown by Ref. 7], the concurrence is a continuous func- 
tion of the anisotropic parameter and reaches its maxi- 
mum at the transition point. Therefore, in both cases, 
one observes that the concurrence itself manifests inter- 
esting behaviors at the quantum phase transition points. 
However, we should emphasize that, such behaviors alone 
may not always signal a quantum phase transition, as 
pointed out by the authors of Refs. 0,13- 

In Ref. 1^, we studied extensively nearest-neighbor 
spin entanglement in the antiferromagnetic XXZ chain. 
By applying results derived from the Bethe ansatz so- 
lution of the model, we showed clearly that the concur- 
rence between two spins located on a pair of nearest- 
neighbor sites in the system is a continuous function of 
the anisotropic coupling parameter and becomes maxi- 
mal at the isotropic Heisenberg point. In this paper, we 
continue our discussions on this issue. Our main pur- 
pose is to show that some fundamental properties of the 
XXZ model, such as non-degeneracy and concavity of 
the ground state energy of the system at the phase tran- 
sition point, commands strongly on the behavior of the 
concurrence. Therefore, we expect that the same scenario 
will appear in a wide class of localized spin models, such 
as the spin ladder model and^ in particular, the XXZ 
model in higher dimensions It is well known that, 

as far as the above-mentioned properties are concerned. 



the ground states of these models are akin to the antifer- 
romagnetic XXZ chain. 

This paper contains two parts. In the first part, based 
on some well-known facts about the antiferromagnetic 
spin models, we prove rigorously that, when the anti- 
ferromagnetic XXZ model is defined on a d-dimensional 
finite bipartite lattice, the concurrence between two spins 
located on a pair of nearest-neighbor sites is an analyt- 
ical function of the anisotropic parameter and takes on 
its maximum at the Heisenberg isotropic point. Then, in 
the second part of this paper, justified by the existence 
of magnetic long-range order (LRO) in the XXZ model, 
we use the spin-wave theory to show that a cusp-like be- 
havior of the concurrence develops in the thermodynamic 
limit when the dimensionality of the lattice d > 2. 

To begin with, we first introduce several notations. On 
a finite d-dimensional simple cubic lattice A with N\ = 
L'* sites, the Hamiltonian of the antiferromagnetic XXZ 
model is 
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where Sf, S^ and S? are spin-1/2 operators at site i 
and A = Jz/ Jx [Jx = Jy) is a dimensionless parameter 
characterizing the anisotropy of the model. The sum in 
the Hamiltonian is over all pairs of nearest-neighbor sites 
i and j. Obviously, this Hamiltonian commutes with the 
total spin z-component operator S^^^^^ = 

Ei^f- Thus, 

each eigenstate of the Hamiltonian is also an eigenstate 
of S^^^^i- Consequently, the Hilbert space of the system 
can be decomposed into numerous subspaces V{M). In 
each subspace, the spin number S^^^^^ = M is specified. 
It is well known that, on a finite simple cubic lattice A, 
the ground state of the XXZ model is nondegenerate in 
any admissible subspace V(M) [isl IT^. In particular, 
its global ground state 5*0 (A, A), which coincides with 
the ground state of the model in the subspace V{M = 0) 
fl^, is also nondegenerate. Therefore, all the physical 
quantities, such as the ground state energy -Bo(A, A) 
and the spin correlation function {S?S?) are analytical 
functions of the parameter A, as long as the lattice is 
finite. 

The conservation of i^totai implies also that, with re- 
spect to the standard basis vectors | IT), | 11), \ it) and 
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I II), the reduced density matrix of two spins on a pair 
of nearest-neighbor lattice sites i and j can be put into 
the following block-diagonal form 



Pij = 




(2) 



As a result, the concurrence of the two spins is Cy ~ 
2 max l^l^l — Vu+u^, 0^ terms of the correlation 

function Gf" = {S"S"), a = x,y,z, it can be exphcitly 
written asfisj 



2 max 
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By the variational principle, one can show that all the 
spin correlations functions G^" are negative. Thus, one 
has 
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-e?j(A, A) - i + (A - l)Gr/ 
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where e?j(A, A) = Eq{K, A)/Nb {Nb is the number of 
bonds in the lattice) is the ground state energy density 
per bond. Furthermore, since all quantities in Cy are 
analytical functions of the parameter A, we are allowed 
to take derivatives of it with respect to A. In particular, 
after taking the first order derivative of Cy , we obtain 



dA 



dA 



+ (A - 1) 



dA 
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Again, due to the nondegeneracy of the global ground 
state \E'o(A, A) of the XX Z model on a finite lattice, we 
can use the Hellman-Feynman theorem to calculate the 
derivative 9e?j(A, A)/9A, which equals to Gfj^. There- 
fore, we finally obtain 



dA 



2 (A 



1) 'J 



2(A-1) 



OA 
g^4(A, A) 
9A2 



(6) 



Immediately, one sees that A = 1 is an extreme point of 
the concurrence. 

Next, we show that A = 1 is actually a maximal point 
of Gij and the concurrence does not have other extreme 
point. In fact, both the statements are the corollaries 
of concavity of the ground state energy Eo{A, A) of the 
Hamiltonian Hxxz with respect to the anisotropic pa- 
rameter A. By the variational principle we know 
that, for any two parameters Ai and A2, the following 
inequality 

E„iA, AAi-f (1-A)A2) > \Eo{A, Ai) + (l-A)So(A, A2), 

(7) 



where < A < 1, holds true for the ground state energy 
-Eo(A, A). In particular, when Eq{A, A) is differentiable 
with respect to A, inequality is equivalent to 



d^EpjA, A) 
dA^ 



< 0. 
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Consequently, we have also (9^6^ (A, A)/dA'^ < 0. Now, 
let us take derivative of Eq. © again with respect to A. 
It yields 



9A2 
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Therefore, A = 1 is indeed a maximal point of the con- 
currence. 

Finally, we prove that A = 1 is the unique extreme 
point of the concurrence Gy. For that purpose, we no- 
tice that inequality (|HJl is actually strict. In other words, 
the equal sign in it can be ignored. This can be easily un- 
derstood by observing the following fact: As A increases 
from -00 to 00, quantity {S[S[) = de9-{A, A)/dA be- 
comes more and more negative. Consequently, the prod- 
uct on the right-hand side of Eq. (O cannot be zero at 
any point except A = 1. That completes our discussion 
on the general behavior of the concurrence Gy for the 
antiferromagnetic XX Z model on a finite d-dimensional 
simple cubic lattice. In addition, we point out that the 
above proof can be easily extended to other cases, such 
as the spin ladder model at J± = 0. 

In the following, we shall discuss the behavior of the 
concurrence in the thermodynamic limit. In Ref. 0, by 
using the Bethe ansatz solution of the one-dimensional 
XXZ chain, we obtained the explicit expression of the 
concurrence near the isotropic point 



Gm+1 = Go-Gi(A-l) 



(10) 



where Go and Gi are two real constants. Therefore, the 
concurrence of the ID XXZ chain is a differentiable func- 
tion of A in the thermodynamic limit. However, things 
are quite different in higher dimensions. For the XXZ 
model in higher dimensions, there exists no exact solu- 
tion. One either uses approximate analytical approach 
such as the spin- wave theory or numerical approach such 
as exact diagonalization studies of finite lattice. To ob- 
tain results in the thermodynamic limit, finite size scal- 
ing analysis must be performed. By using the stochastic 
series expansion quantum Monte Carlo method for lat- 
tices up to 16 X 16, Sandvik '2^ did an extensive study 
on the two-dimensional S=l/2 antiferromagnetic Heisen- 
berg model. The finite size results for various ground 
state quantities were extrapolated to the thermodynamic 
limit using fits to polynomials in 1/L, constrained by 
scaling forms previously obtained from renormalization- 
group calculations for the nonlinear sigma model and 
chiral perturbation theory. He demonstrated that the 
results were fully consistent with the predicted leading 
finite size corrections. With the same scaling forms, Lin, 



Flynn, and Betts [21| studied the XX Z model on square 
lattices and obtained various quantities as functions of 
anisotropic parameter A for the infinite system. Two 
conclusions from pervious works mil 122", "2^ are rele- 
vant to the present study: (i) results obtained from the 
spin-wave theory are qualitatively correct and quite ac- 
curate, usually within 3 percent as compared with exact 
solution on finite lattices; (ii) derivatives of the ground 
state energy with respect to the anisotropic parameter 
A are not continuous at the Heisenberg point A = 1, 
for example, see Figure 3 in Ref. This conclusion 

is consistent with the belief that there exists antiferro- 
magnetic long-range order (LRO) in the d-dimensional 
XXZ model for d > 2. In fact, the existence of the 
LRO for d > 3 was rigorous proved ^24*1, while for d = 2 
most numerical studies support it. Based on these con- 
clusions, we apply the spin-wave theory to calculate the 
concurrence Cy of the XXZ model. We also use exact 
diagonalization results as complementary. As shown in 
the following, the symmetry breaking in the thermody- 
namic limit, which is absent in the one-dimensional case, 
causes the singular behavior of the concurrence at the 
quantum phase transition point. 

Following the standard procedure, the XXZ Hamil- 
tonian is mapped into a boson model via the Holstein- 
Primakoff (HP) transformation 



5*;+ = V2S{l-hi/2S) 
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'^{l-fii/AS) di, 



Sr = V2Sdl{l - ni/2Sf'^ ~ V2Sdl{l - fn/AS), 
S? = S-dldi, (11) 

where di and a[ are boson creation and annihilation oper- 
ators at site i for the spin deviation. In the region A > 1, 
the antiferromagnetic ordering is in the spin-z direction. 
Consequently, we have 

Hxxz/^ ^ ^ [S^ + S (dldi + ajflj) 



(ij> 

-|- xS ( didj + a[a| 



(12) 



where x = 1/A. Using Fourier transform, we rewrite the 
Hamiltonian as 

Hxxz/A = -^NS^ + zSj2H{-k), (13) 



where z is the coordination number of the lattice and 

(14) 



H(k) ^ aj^flk + (^akft-k + a_k"k 



with 7k = I J2m=i coskm- By applying the Bogoliubov 
transformation 



Sk = WkCk - WkcLk 

Ok = -WkC-k + WkCk, 



0.25 



0.15' 



0.05 



2D XXZ model 




B-^ 4*4 
6*6 

— Spin wave theory 



FIG. 1: The concurrence Cij of the two-dimensional XXZ 
model as a function of A(= Jz/Jx)- In the Figure, the 
doted lines are obtained from the exact diagonalization for 
4 * 4(square) and 6 * 6(circle) square lattices respectively. 



we diagonalize -ff (k) and obtain 

i/(k) = Wk - 2;7kUkWk 

+ (uk + ^'k - 2xukfk7k) cj^Ck, (16) 

where the Uk and Vk satisfy the following constraint con- 
ditions 



2 2 1 

a;7k. 2 , 2n n 

-— (Uk + I'k) - "kWk = 0. 



(17) 



Finally, the ground state energy of the model in the re- 
gion of A > 1 can be written as 

EoiA > 1) = -'-NS' + ^ l) -(18) 

By similar approach, we can also obtain the ground 
state energy of the XXZ model in the parameter region 
of < A < 1. In this case, the system has antiferro- 
magnetc order in the XY plane in the thermodynamic 
limit. As a result, the diagonalized Hamiltonian has the 
following form 

iJ(k) = (1 + y7k) vl ~ x7kUkWk 

+ [(1 -I- y-fk) (uk + Wk) - 2a;ukWk7k] 4ck(19) 

where x = (1 + A)/2 and y — {1 — A)/2, and the corre- 
sponding ground state energy is 



En 



zS ■ 
~2 



iV52 + ^^(l + y7k) 



1- ^'7^/(1 + 2/7k)2-l) 



(20) 



Within the spin-wave theory framework, we calculate 
the spin correlation function G\? and hence the concur- 
rence Cij of the model in two and three dimensions. Our 
(15) results are shown in Figs.^andl^l respectively. We also 



4 




FIG. 2: The concurrence Cij of the three-dimensional XX Z 



model as a function of A(= Jz/Jx)- 



show results obtained from the exact diagonahzation of 
the XXZ model on finite square lattices. The trend as 
function of lattice size is clear. It is interesting to see 
that, in both cases, the concurrences Cy of the XXZ 
model not only have their maximal value at the critical 
point A = 1, but also show discontinuities in their first 
derivative with respect to A at the transition point. This 
behavior is quite different from the one-dimensional case 
(Eq. I10|l . as we expected. We attribute this difference 
to the existence of the magnetic long-range orders in the 
system with d> 2. 

As we have seen, the concurrence Cy is closely related 
to the ground state energy of the model. As a result, any 
singularity in the ground state energy may be inherited 
by the concurrence On the other hand, on a finite 
d-dimensional simple cubic lattice, the ground state of 
the antiferromagnetic XXZ model is non-degenerate for 
A G (— oo, oo),16j. Therefore, the ground state energy 
Eq{A, A) as well as the concurrence Cy are analytical 
functions of A, regardless of the dimensionality of the lat- 
tice. However, it is no longer true in the thermodynamic 
limit. For the one-dimensional XXZ model, it is well 



known that its ground state in both A < 1 and A > 1 
regions does not have magnetic long-range order. There- 
fore, we do not expect a dramatic change in the ground 
state energy Eo taking place at A = 1. Consequently, 
the concurrence will behave more or less like itself on a 
finite lattice. However, in two and three dimensions, the 
ground state energy of the system develops a cusp at the 
transition point in the thermodynamic limit 21]. This 
phenomenon can be understood by the picture of the 
first excited-energy levels crossing at A = 1, required by 
the existence of magnetic long-range order [2^. There- 
fore, a singularity inherited by the concurrence at the 
transition point, is expected to appear, as shown by our 
calculations. 

In summary, we have studied the ground state two-spin 
entanglement, as measured by the concurrence, in the d- 
dimensional XXZ model. We gave a rigorous proof that 
the ground state concurrence in the XXZ model reaches 
maximum at the isotropic point. We extended our previ- 
ous studies in one dimension j3| to two and three dimen- 
sions by using the spin-wave theory and exact diagonal- 
ization technique. The use of the spin-wave theory is jus- 
tified by the existence of magnetic long-range order in the 
XXZ model for dimensionality d> 2. We found that the 
concurrence in two- and three-dimensional XXZ mod- 
els also reaches maximum at the isotropic point A = 1. 
Unlike the one dimension case, the concurrence shows 
cusp-like behavior around the critical point, and its first 
derivative is not continuous in the vicinity of the critical 
point. 
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